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lit is known that arrays of trapped ions can be used to efficiently simulate a variety of many- 
body quantum systems. Here, we show how it is possible to build a model representing a spin 
chain interacting with bosons that is exactly solvable. The exact spectrum of the model at zero 
temperature and the ground-state properties are studied. We show that a quantum phase transition 
occurs when the coupling between spins and bosons reaches a critical value, which corresponds to 
a level crossing in the energy spectrum. Once the critical point is reached, the number of bosonic 
excitations in the ground state, which can be assumed as an order parameter, starts to be different 
from zero. The population of the bosonic mode is accompanied by a macroscopic magnetization of 
the spins. This double effect could represent a useful resource for the phase transition detection 
since a measure of the phonon can give information about the phase of the spin system. A finite- 
temperature phase diagram is also given in the adiabatic regime. 
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I. INTRODUCTION 

Quantum spin chains play a fundamental role in the 
study of many-body systems and quantum phase transi- 
tions. These phenomena take place at zero temperature, 
that is in a purely quantum regime, and are induced by 
the variation of an internal parameter causing a critical 
change in the ground state due to level crossing in the 
energy spectrum In recent years, a renewed interest 
in quantum phase transitions has been developed to un- 
derstand the behavior of entanglement near the critical 
point [2 5] . On the contrary, the study of interactions be- 
tween spins and light is considered as the starting point 
to introduce the framework of open systems, where the 
appearance of dissipation and decoherence can be under- 
stood A prototypical example is the so-called spin- 
boson model 0, consisting of a single spin interacting 
with a multimode electromagnetic radiation, modeled as 
distribution of quantum harmonic oscillators. Here, we 
discuss the case where a series of spins (arranged in an 
isotropic XY chain) interacts with one or more bosonic 
modes. 

Despite the importance of spin systems, they are not 
always directly accessible experimentally, and efficient 
experimental quantum simulation methods are needed to 
study their properties. It has been shown that spin chains 
can be efficiently simulated using internal and external 
degrees of freedom of trapped particles [HlH . In partic- 
ular, arrays of laser-cooled trapped ions seem to be very 
promising from an experimental point of view. Ions can 
be trapped with high spatial accuracy, and their internal 
states can be manipulated with high precision by means 



of the interaction with electromagnetic fields. The first 
experimental evidence of coupling between two-level sys- 
tems, consisting in the transition from paramagnetic to 
ferromagnetic order in a tw o-sp in quantum Ising model 
has recently been reported [l3[. Furthermore, ion traps 
offer the possibility of engineering spin-boson coupling 
[ill Il4j . Since spin chains are also important for practi- 
cal applications, like quantum communication protocols 
[Tlj . their simulation would allow estimation of the de- 
coherence effects which derive from the interaction with 
the environment [l6j]. 

As said before, we propose a model which allows to ob- 
serve the phase transition in an isotropic XY spin chain 
coupled with external boson modes. The model we in- 
troduce in the following is exactly solvable. The phase 
transition manifests itself in a nonanalytical variation of 
the amplitude of the bosonic field. In correspondence 
of this change, the chain acquires a finite magnetization. 
The physical implementation of this model can be done 
within the framework of trapped ions discussed before. 

The paper is organized as follows. In Sec. [H]wc first in- 
troduce the general argument of simulation of spin chains 
with trapped ions. Then, we show how it is possible to 
simulate a particular model, consisting of a chain of spins 
in contact with a bosonic environment. The critical prop- 
erties at zero temperature of this model are discussed in 
Sec. IIII1 while in Sec. |TV]we study the phase diagram at 
finite temperature in the case of a single slowly oscillating 
phonon. Finally, in Sec. |V]we present our conclusions. 



II. ARRAYS OF TRAPPED IONS AND 
INTERACTION BETWEEN SPINS AND BOSONS 
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As shown in Refs. @ and 0], Coulomb chains in linear 
Paul traps can simulate spin-spin interactions. An effec- 
tive spin-spin Hamiltonian emerges as the result of the 
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interaction of the internal (electronic) states of ions with 
the phonemic degrees of freedom generated by Coulomb 
repulsion. The total Hamiltonian includes a phonon bath 
(H u ), a state-dependent force produced by a set of lasers 
along the directions a = x,y,z (Hf), and an effective 
magnetic field that can be generated by forcing transi- 
tions between the internal ion states (H m ). In units of 
% = 1, we have (see @) 



H f = -2£)F a?a> ,(l + <7f), 

a,l 



we obtain 



(1) 



where a Qi „ is the annihilation operator of the n vibra- 
tional mode in the a direction of the phonon bath due to 
the Coulomb repulsion, F a are the laser forces, q a j is the 
displacement, with respect to its equilibrium position, of 
ion I in direction a, and B a is the effective magnetic field. 
Expressing the coordinates in terms of collective modes 
through the matrices M a , 



a,l.n 



y / 2m6Ja^ 



«n + a«,«)(l + of). (2) 



A spin-spin Hamiltonian is obtained after the applica- 
tion of a suitable canonical transformation U — e s , in- 
troduced to eliminate Hf-. if 



a,Ln 



= «n-aa.«)(l + of), (3) 



the total system Hamiltonian turns out to be 

H = UHU^ = H„+ l -Y, J U< a f + E B ' a < + H ^ 

a>.l,j a, I 

(4) 

where B' a = B a — ^/(muj), the coupling parameters 
J are function of F, M, ui, and He is a residual term that 
can be neglected at low temperatures or by using highly 
anisotropic traps @]. A particular case of Eq. (|4]) is 
represented by the isotropic XY chain in the presence of 
a transverse field, numerically studied in Ref. @ . 

Here, we show that a different interpretation of the 
unitary transformation allows us to obtain other kinds of 
phase transitions, involving, for instance, the interaction 
between spins and bosons. A possible way to do that 
consists in applying different transformations along the 
three axes. In fact, by limiting the transformation given 
before to the directions x and y, through 



S x ,y 



a=x,y l,n 



M, 



I , n 



= (oU-Oa,»)(l + *■?)» (5) 



M? 



= (' I L+ fl z,n)( 1 + ff l) 



l,n V 2muJ z,n 

+\ E Ew+E^ 

a=x,y IJ I 

+ E Y, B ' a(7 " + H E+ H ^ ( 6 ) 

a—x^y I 

Here, the correction H' E can be neglected in the same 
limit of He- In the case of a spatially homogeneous array, 
we can drop the dependence from the ion position in both 
M and in J. Moreover, if the Coulomb interaction can be 
considered as a perturbation with respect to the trapping 
potential (stiff limit), the trapping frequencies can be 
tuned to obtain relevant spin-spin coupling only between 
nearest neighbors Q. Along the transverse direction, a 
boson displacement can be applied to eliminate the term 
E„ F zMZ/(y/2mu} z , n )(at, n + a z . n ). As a result, 

6 = E ?E«i+E^««+ a ^)E^ 



a—x,y 



z2 



+(*"-2£ff-)E°f + E E B ' a(j i 

n 71 I a—x,y I 

— * /,}.. 



(7) 



where = F z M*^y '(N /2mu Zi „). The last term is a 
constant that translates all the energies, and the fields 
B' can be set to zero. By fixing the external magnetic 
field B z = 2^2 n (g^ /w n ), we arrive to 



iV 



a z,n "r a z,n 



- NB Z , N 

)£of + (8) 



It is worth noting that, since the mutual distance between 
the ions makes any direct spin-spin interaction negligible, 
the spins we introduced are effective two-level systems 
that, because of the unitary transformation, are not triv- 
ially related to the original internal degrees of freedom 
of the ions. Nevertheless, the bosonic variables do re- 
fer to the real phonon modes along the z direction, so 
their measurement appears to be the natural way to get 
information about the overall system. 



III. CRITICAL PROPERTIES 

To describe the critical properties of this model, let us 
first discuss the simpler single-mode version. By assum- 
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ing J x = J y 

N 



J, the Hamiltonian is 



N 



H sm = J^2lrf<7 l+1 +H.c.]+uja i a+-j=(a' ! + a)^Vf, 
1=1 J=1 

. (9) 

where N is the total number of spin, and where periodic 
boundary conditions are imposed. As it is useful to work 
with adimensional quantities, we use J as energy unit 
and define the rescaled phonon frequency 7 = lj / J as 
well as the bare coupling constant A = g 2 /(Jui). Using 
this notation the Hamiltonian reads 



N 

H sm = YJ a i~ a r+ 1 + H - c -\ + 7 flta - 

1=1 



N 



-0 + + a)]>>f 



1=1 



(10) 

To diagonalize it, first we introduce the Jordan- Wigner 
transformation (it| . defined by erf = 1 — 2c|q, cr ( + = 

I1 J<; (l - 2ctcj) q, and crj" = IIj<i (l - 2c ] c j) 4> 
mapping spins into spinlcss fermions. Then, the Fourier 
transform 



Ck 



N 

E- 



(11) 



i=l 



allows us to write, in the thermodynamic limit, 




2c\c k ) 



H sm = Y e kcl°k + 1°^ a + 
k 

(12) 

with Cfe = — 2 cos(27rfc/iV). Fermions can be decou- 
pled from the bosonic degree of freedom by applying 
the displacement operator D = cxp [a{a^ — a)} such that 

DaD^ = b,a= - ^/X/iW) E fc (l - 2c { c k) being an op- 
erator in the fermionic space. Since [d, ^2 k ekC k Ck] = 0, 
fermions are not transformed into polarons because of the 
application of D. Moreover, it is immediately seen that 
[b, b^} = 1. This displacement replaces the boson vacuum 
with a coherent state whose amplitude is determined by 
the transverse component of the total spin. The final 
Hamiltonian is 



H, 



E 



tfb 



-Y 

k 



(1 - 2c\c k )] 2 . (13) 



Adding the boson displacement D to S z and S x , v of Sec. 
IIIl amounts to building the total S of Ref. Q for a par- 
ticular choice of parameters. If one disregards the term 
jb^b. the Hamiltonian Q2] can now be interpreted as a 
XX Z spin chain with short-range coupling along the ra- 
dial direction and axial long-range coupling. These kinds 
of interaction ranges are realistic, since radial modes 
can mediate nearest-neighbor spin-spin interactions (stiff 
limit), and transverse modes are known to generate long- 
range interactions (soft limit). While the range of in- 
teractions along the transverse direction depends only 
on the number of ions and cannot be manipulated by 



means of laser fields, the stiff limit in the radial direction 
can be reached by increasing the trapping frequencies 0] . 
Performing the unitary transformations in two separate 
steps helped us to interpret this model in terms of the 
coupling with a real phonon. Without performing the 
two steps separately, we could not discuss the spin-boson 
phase transition. 

At first sight, H^ m is similar to that of the XX chain in 
a transverse field [18( , whose Hamiltonian in the Jordan- 

Wigner space would be J2k ekC l Ck ~ ^Sfc(l — 2c \ c k)- 
While the latter model has the U(\) symmetry, since it is 
invariant under rotations around the z axis, the Hamilto- 
nian ^ has an additional symmetry, also being invariant 
under the action of the operator S = Yli of <S>exp [«7roTa]. 
The U(l) symmetry implies, in the language of fermions, 
that eigenstates of H have a fixed number of particles. 
The symmetry S could be broken, for instance, by adding 
a displacing bosonic field, a + a} , and then sending it to 
0. The average value of a + can be assumed as an 
order parameter to study the breakdown of S. 

Since we are interested in the implementation in an ar- 
ray of ion traps, finite-size effects should be taken into ac- 
count. A detailed analysis of XYZ models in finite cycles 
has been made in several papers (l8l - [20l |. In our model, 
after the Jordan- Wigner transformation, the Hamilto- 
nian becomes 

N-l 

H sm = Y, ( c '+i4 + c ' c l+i) ~ ^{ c i c n + cjyCi) + 7a t a 
i=i 

+y^(« t +«)E( i - 2c ^), (14) 

k 

where N is assumed to be an even number and where 
V = rL=i(l — 2cJq). Its possible eigenvalues are ±1. 
Note that V is a measure of the parity of the number 
of particles of each state. Since [H sm ,V] = 0, all eigen- 
states of H sm have definite parity, and we can proceed to 
a separate diagonalization of H sm in the two subspaces 
corresponding to V = ±1. Then, if we introduce the 
two new Hamiltonians H^ m = H sm (V = ±1), the com- 
plete set of eigenvectors of H sm will be given by the odd 
eigenstates of H~ m and the even eigenstates of H+ m . The 
choice oiV = — lor7 :, = +l amounts to having, respec- 
tively, periodic or antipcriodic boundary conditions. This 
implies that Hf m and H~ m are diagonalized through two 
Fourier transforms which differ from each other because 
of the set of allowed values of k. If half- integer values 
(k = 1/2, 3/2, . . . , N — 1/2) are used to diagonalize Hf m , 
in the case of H~ m we have k = 1, 2, . . . , N. Now, the 
displacement operator D can be introduced, and a struc- 
ture formally identical to that given in Eq. (|13|) can be 
done for both H+ m and H~ m , with the proper choice of 
k. 

The ground state of the bosonic part can be identi- 
fied as the vacuum state of the operator b, and it corre- 
sponds, in the original representation, to a coherent state 
of amplitude equal to the average value of the operator 
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a calculated on the fermionic ground state, whose struc- 
ture is now enlightened. Since both Hf m and H~ m are 
exactly solvable, their spectra can be calculated by mea- 
suring the energy of all possible configurations, obtained 
adding the desired number of fermions. The energy of a 
state with m particles is 



-5 - 



w -10 



Em — - 



X(N-2m) 2 2nk 

N - 2 E C0S — ■ (15) 



where M. is a string representing the occupation of differ- 
ent modes. The sum is extended only to the values of k 
where a particle is present. For large values of A, the first 
term will be larger than the sum, unless TV — 2m is very 
small. The ground state corresponds to the sequences 
which maximize (TV — 2m) 2 , that is, to m = or m = TV, 
and it is twofold degenerate. In fact, both the state 
|$+) = |0) ® | a) and the state |$") = n fc 4l°) ® I ~ «) 
have energy equal to Eq = —NX. Here the state |<& + ) is, 
in the original fermion-boson representation (before the 
application of D), the tensor product of the fermionic 
vacuum and of a coherent bosonic state of amplitude 
a = —gy/N/oj, while in |$ _ ) the fermionic system is 
fully occupied, and the boson coherent state has ampli- 
tude —a. As A decreases, the two terms in Eq. (fl~5j) start 
to compete with each other. Let us call A m the value such 
that Eq = min{i?x}. Since, by symmetry, A m = Aj\r_ TO , 
we limit out consideration m = 1,2,..., N/2. It can be 
shown that, if m < m' , then A m < \ m >. In other words, 
the first ground-state level crossing takes place between 
|4>+) (or |$~)) and a state with N/2 particles. This tran- 



sition occurs at the value of A 



N/2 



A c determined by 



AT/4 



TV A, 



cos 



fc=l 



(2k- l)?r 
TV ' 



(16) 



In the thermodynamic limit, performed replacing the 
sum with an integral, we obtain A c = 2/tt. As a con- 
sequence of the half-filling, ^2 k (l — 2c^Cfe) = 0, and the 
boson part is left in its vacuum. This state can be written 

as 1*^) =111^/4 4l0>®|0)- 

By further decreasing A, no other transitions are ob- 
served. Then the state \$> HF ) is the Hamiltonian ground 
state for every A < A c . It is worth noting that both 
\$ HF ) and |$ + ) (or |$~)) are eigenstates of H+ m . In 
any case, eigenstates of H~ m are excitations. In Fig. [TJ 
we plot the lowest energy levels for a chain of eight spins. 
The transition — > \& HF ) is observed. 

While \$> HF ) is an eigenstate of the symmetry opera- 
tor 5, this is not true in the case of |$ + ) and |$~). In 
fact, 5|<I> ) = The possibility of obtaining a coher- 

ent emission of light on the bosonic mode is inherently 
linked to the breakdown of the Hamiltonian symmetry, 
since the only coherent state which is also an eigenstate 
of exp [ina^a] is the vacuum. While symmetry breaking 
takes place independently of the system size (the num- 
ber of spins), only in the thermodynamic limit docs a 
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FIG. 1: Ground-state energy of H as a function of A for a 
chain of 8 spins. Both the energy and A are measured in units 
of J or a; (we assume J = ui = 1). For any given number of 
particles, only the lower level is plotted. The solid green line 
represents the energy of |$ + ) and |$~), while the solid black 
line corresponds to \<b HF ). As we can see, there is only one 
transition, occurring at A c ~ 0.65. All energies derived from 
numbers of fermions different from 0, TV/2, TV are excitations. 
Specifically, the dotted green line represents the energy of the 
one-particle state, the dotted red line is for the two-particle 
state, and the dotted black line corresponds to the three- 
particle state. Due to the particle-hole symmetry, states with 
m and TV — m particles have the same energy. 



true phase transition take place, since the Hilbert spaces 
that can be built up starting from |<f> + ) and |$~) become 
unitarily nonequivalent. 

The exact results we propose have been obtained in the 
presence of periodic boundary conditions. By releasing 
this hypothesis, they are correct only in the thermody- 
namic limit, while in the case of finite size, as in the 
example in Fig. [TJ corrections to the energy levels are 
expected. Actually, the structure of eigenfrequencies giv- 
ing rise to the phase transition from \^) to \Q HF ) can 
be observed also in the case of a very short open chain 
(TV = 4). 

In the presence of multimodc radiation, Hamiltonian 
(J5J) can be recovered modifying Hamiltonian © in the 
following way 



TV 



1 = 1 



N 



-EV^^ + ^E^, (17) 



i=i 



and it is subject to the symmetry S n 



exp[i7rX; n 4, 



The diagonalization is performed 



through the application of the Jordan- Wigner trans- 
formation and of the operator FJ D n , with D n = 

exp[d n (al-a n )], and d n = - V / A„/(7„TV) £) fc (l - 



5 



2ctc k ) 



As a result, we have 



A 



H mm = E £kc{c k + E Unb^bn - JjlJ2( l - 2c l c k)} 2 , 
k n k 

(18) 

where A = J2 n 9nl '(J^n) is related to the spectral den- 
sity of the bath. In Ref. [IJ], the authors showed how 
baths characterized by different spectral densities can be 
simulated. Being H mm formally identical to H sm , we ex- 
pect a phase transition for A c = 2/tt. For A > A c , the 
two degenerate ground states are |$+ m ) = |0) ® T7 n l a «) 

and |$ mm ) = n*4l°> ® Il« I - «n>. while for A < A c 
the non degenerate half filled ground state is = 

n\k\<N/A\°)®u n \o)- 

To clarify the proposal, we comment on the prepara- 
tion scheme needed to observe phase transitions in this 
system. As mentioned in Q, the procedure consists in 
initialization to the state |1 ... \.}, and then adiabatic 
switching of the spin-spin couplings. This process is re- 
peated several times for different ratios of J/g, while the 
phase transition is detected through fluorescence of the 
ions internal levels, as experimentally realized in [l3j |. 
Here, we also need to perform motional ground-state 
cooling as far as axial degrees of freedom are concerned. 
The possible errors in simulating this model with trapped 
ions comes from the canonical transformation and H' E , 
where the error is quantified by the boson displacements 
and the mean phonon number (given by temperature) of 
the neglected transverse modes. Therefore they have to 
be cooled down too I8f . 



IV. TEMPERATURE EFFECTS IN THE 
ADIABATIC LIMIT 



In this section we analyze the phase transition at finite 
temperature in the so-called adiabatic limit (7 1) with 
a single coupled mode. This limit corresponds to the case 
of a very slowly oscillating phonon and it is asymptoti- 
cally exact at high temperatures. The adiabatic Hamil- 
tonian is obtained by neglecting the kinetic energy of the 
phonon and treating the coordinate as a parameter 



N 



H AD =J2i (J t a i+i+ H - c - 



1=1 



v 

4A 



N 



E CT ?' (19) 



(=1 



where we introduced the adimensional coordinate v = 
(g/ J)\/2uimx. The Hamiltonian is a isotropic XY model 
in an external field, parametrically dependent on v, it can 
be diagonalized in the same way as described previously 
and reduced in two blocks corresponding to odd and even 
pscudofcrmion occupation number 



#4 



E^ 



{v)c\c k 



V 

4A 



N 



(20) 




FIG. 2: Phase diagram in the adiabatic approximation. Data 
were obtained with 10 spins and numerical evaluation of the 
adiabatic potential. The left region corresponds to a single- 
well adiabatic potential, the right region to a double-well po- 
tential. 



(v) = -2 cos(2irk/N) +v/VN 



with 



of a configuration with a filling m is 



The energy 



E 



M 



(vm \ ^ 2ir \ 
—= + > cos — k 
v keM / 



(21) 



The partition function, defined by 



dvZ {v,P) exp-/3(^ + -)=), 



(22) 

with the adimensional parameter ft = J/ksT and where 
ZoM=tr{e-^* D ^}, (23) 



can be written as Z oc J duexp — ftVAD^v, ft) so defining 
an adiabatic potential 

V AD (v,ft) = ^ + -?= + -~lnZo(v,ft). (24) 

Above the critical temperature, Vad iy, ft) is minimized 
by v = 0, while, at the critical point, the potential admits 
two minima, allowing for a displacement of the phonon. 
The phase diagram obtained is plotted in Fig. [2] 

The experimental procedure to measure the phase 
transition is exactly as in the zero-temperature case. The 
radial vibrations should be cooled, while the axial degrees 
of freedom can be prepared in a "thermal state". The 
temperature of this state can be controlled by reducing 
the power of the laser (resolved sideband) cooling mech- 
anism. To raise it further, the power of the Doppler laser 
cooling has to be tuned. We still expect to observe an 
abrupt change in fluorescence properties, as in the zero- 
temperature case, since for moderate temperatures the 
coherent state |a) becomes thermal, but the mean num- 
ber of phonons does not disappear abruptly. According 
to Fig. O the order parameter value should change as a 
function of T . 
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V. CONCLUSIONS 

To resume, we have shown that a system of trapped 
ions can be mapped into an isotropic XY chain inter- 
acting with phonons, following the scheme in Ref. Q. 
Modifying the original canonical transformations needed 
to write a spin-spin Hamiltonian, the vibrational degrees 
of freedom along a fixed direction are coupled with the 
effective spins the other degrees of freedom are mapped 
onto. The resulting model is exactly solvable and exhibits 
a quantum phase transition, due to the breaking of a sym- 
metry which takes into account both the spin and the bo- 
son degrees of freedom. In this phase, phonons are in a 
coherent state with finite amplitude. Since the phonons 
refer to the real ionic vibration and not to an effective 



quantity, phase transition detection should be possible 
using state-of-the-art techniques. The experimental re- 
alization of the building block of such simulations (Ref. 
[IH), encourages this kind of investigation. 
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